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. Abstract 

X 

^ ' We discuss the two- and three-point correlators in the two-dimensional three-state Potts model 

in the high-temperature phase of the model. By using the form factor approach and perturbed 
conformal field theory methods we are able to describe both the large distance and the short 
distance behaviours of the correlators. We compare our predictions with a set of high precision 
Monte-Carlo simulations (performed on the triangular lattice realization of the model) finding a 
complete agreement in both regimes. In particular we use the two-point correlators to fix the 
various non-universal constants involved in the comparison (whose determination is one of the 
results of our analysis) and then use these constants to compare numerical results and theoretical 
predictions for the three-point correlator with no free parameter. Our results can be used to 
shed some light on the behaviour of the three-quark correlator in the confining phase of the 
(2-|-l)-dimensional SU(3) lattice gauge theory which is related by dimensional reduction to the 
three-spin correlator in the high-temperature phase of the three-state Potts model. The picture 
which emerges is that of a smooth crossover between a A type law at short distances and a Y 
type law at large distances. 



1 Introduction 



The aim of this paper is to study the three-point correlator in the Z3 Potts model outside the critical 
point. In particular we shall study the thermal perturbation of the model, in the high-temperature 
phase (i.e. the phase in which the symmetry is unbroken and the magnetization is zero). We shall first 
discuss the two-point function and compare it with high precision numerical simulations in order to 
fix all the normalizations and the non-universal quantities which appear in the three-point function. 
Then wc shall address the three-point function and compare it again with numerical simulations. 
Thanks to the preliminary analysis of the two-point function this comparison will not require any 
fitting procedure but will be a direct and absolute comparison between theoretical predictions and 
numerical results. 

On the theoretical side we shall study both the two- and the three-point correlators with two 
different tools. 

• The form factor approach which is essentially a large distance expansion. This approach requires 
exact integrability, a condition which is indeed fulfilled by the scaling Z3 Potts model. 

• The perturbative expansion around the conformal fixed point which is essentially a short dis- 
tance expansion and is completely general, meaning that no specific integrability property of 
the perturbation under study is needed. 

Prom a theoretical point of view this is a rather interesting challenge: 

• In the large distance regime (where the form factor approach is expected to hold) the strong 
coupling expansion suggests the existence of an additional midpoint, inside the triangle spanned 
by the three spins of the correlators, where the strong coupling paths emerging from the three 
spins converge and join. The appearence of this new point is a novelty with respect to the well 
known form factor calculation for the two point functions. The way in which it is obtained is 
non-trivial and is one of the interesting features of our analysis. 

• Similarly it is the first time that the approach of the perturbative, infrared safe, expansion 
around the conformal solution is used for a three point function. This required some non- 
trivial extension of the techniques used in previous works on the two-point functions. 

Prom a physical point of view the scenario which emerges (which strongly resembles what one finds 
when looking at the three-quark potential in lattice gauge theories (LGT) , a correspondence which we 
shall discuss in detail below) is a smooth crossover, as the distance among the three points increases, 
from a short distance behaviour in which the three point function is dominated by the three spin-spin 
interactions along the edges of the triangle to a large distance behaviour in which the strong coupling 
expectation (the three spins joined by a path of minimal length) is fully realized. 

This scenario is confirmed by the numerical simulations which turn out to be in remarkable 
agreement with our theoretical results. Indeed the important consequence of having exact analytic 
results for the two expansions is that we are able to compare our predictions with triangles of any 
size, both smaller and larger than the correlation length. Moreover, even if in the following we 
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shall mainly perform our comparisons in the equilateral case, we can obtain analytic predictions for 
triangles of any shape, thus allowing a very selective test of our results. 

Besides a better understanding of the three-state Potts model a second important motivation 
which we had in mind addressing the present problem is to use our results to shed some light on 
the behaviour of the three-quark potential in SU(3) LGT's. Indeed by dimensional reduction the 
(2+l)-dimensional SU(3) LGT can be mapped into the two-dimensional three-state Potts model and 
in particular the three-quark free energy is mapped into the three-point function of the Potts model 
(we shall discuss in detail this correspondence in sect|7| below). An important open problem in the 
study of the three-quark potential is to understand if the three-quark correlator follows the so called 
"Y" or "A" law (for a discussion of these laws see again sectdbelow). Our results on the three-state 
Potts model suggest that the right picture is a smooth crossover between the two laws. At short 
distance, i.e. for interquark distances smaller than the correlation length, the three-quark potential 
is well described by the A law, which is indeed exact at the critical point (i.e. when the correlation 
length goes to infinity), while at large distances the correct description is the Y law, which becomes 
exact in the strong coupling limit, i.e. for interquark distances much larger than the correlation 
length. This mixed behaviour agrees with the results of some recent simulations performed directly 
in the gauge model [42]. 

This paper is organized as follows. Sect O contains a general discussion of the three-state Potts 
model both on the lattice and in the continuum and its CFT description at the critical point. In 
sect El we discuss the form factor analysis of the large distance expansion both for the two- and the 
three-point functions. Then in sect|3]we shall compare these predictions with a set of high precision 
Monte-Carlo simulations performed on a triangular lattice. SectlHlis then devoted to the study of the 
short distance perturbative expansion which is then compared with the Monte-Carlo data in sectEl 
In sect d we shall then briefiy comment on the implications of our results for the study of barionic 
states in LGT. Two appendices conclude the paper. The first contains some technical steps we omit 
in the main text, while the second is devoted to universal amplitude ratios. 

2 The three-state Potts model 

The two-dimensional, isotropic, nearest-neighbor, three-state Potts model at temperature T on a 
lattice A is defined by the partition function 

{s„} 

with the Hamiltonian 

n = jY.{i-5s^s^,) (2) 

{nn') 

where s„ = 0, 1, 2 are Zs-valued variables on each site^, n G A, /3 = (kBT)^^, and {nn') denotes pairs 
of nearest-neighbor sites [1]. The symmetry group of the Potts Hamiltonian is S3, the permutation 
^For this reason we shall often denote in the following the model as the Z3 Potts model. 
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group of 3 objects. In the following we shall study in particular the case in which A is a regular 
triangular lattice with honeycomb boundary conditions. 

This choice is motivated by the fact that in the following we shall mainly be interested in the 
three-point function of the model, for which a more symmetric choice of arguments is possible on a 
triangular lattice. Notice however that most of our theoretical results are obtained in the continuum 
limit theory and hence, once the non-universal, lattice dependent, normalizations are fixed, hold for 
any lattice A. A general introduction to the Potts model can be found in the review by Wu [2]. 
Recently a set of interesting results (series expansions and exact free energy calculation on strips) 
for the triangular lattice realization of the model in which we are interested appeared in a series of 
papers (see in particular [3], [4] and references therein) to which we refer the interested reader. 

The model is known to have a second order phase transition for a critical value of the tem- 
perature which separates the high-temperature, symmetric phase from the low-temperature one in 
which the symmetry is broken and a spontaneous magnetization appears. In the following we shall be 
interested in the symmetric high-temperature phase. This is due to the fact that (as we shall discuss 
in the last section) we plan to use our results to better understand the behaviour of the baryonic 
states in LGT, and the confining phase of the SU(3) LGT is mapped by dimensional reduction to 
the high-temperature phase of the 3-state Potts model. 

Assuming that a single phase transition point exists in the model, the critical temperature can be 
obtained exactly even for the triangular lattice by using duality and the star triangle relation which, 
combined together, lead to an algebraic equation for Tc [5]. Defining a = we have 

- 3a - 1 = (3) 

whose unique positive solution is 

a = cos(27r/9) + V3sin(27r/9) = 1.879385... (4) 

to which the critical value f3c = 0.6309447... corresponds. 

At the critical point the model is described by a minimal (non-diagonal) CFT with central charge 
c = 4/5. It may be useful to briefly describe this model. It is indeed the simplest example of a non- 
diagonal minimal model of the Virasoro series [6]. However it can be also realized as the simplest 
diagonal minimal model of the so called algebra [7] . 

Its operator content is composed by six primary fields, which however, due to the non-diagonal 
nature of the model, lead to a larger number of operators when the analytic and antianalytic sector 
are combined. Using the standard CFT notation (r, s; r', s') which corresponds to the conformal 
dimensions [hr^sihr' ,s') for the field, with /ir,s given by 

(6r - hsf - 1 

hr,s = (5) 

the most relevant operator in the energy sector is e = (2, 1; 2, 1) with scaling dimension = 2/i2,i = 
4/5, while the most relevant operators in the magnetic sector are the doublet a and a, of type 
(3, 3; 3, 3), with scaling dimension Xa- = X^ = 2/13^3 = 2/15. 

The nice feature of the description is that the operator content is composed by four fields 
only (besides the identity) and all the higher fields appear as secondary fields in the VF3 conformal 
families. 
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3 Large distance expansion 



In this section we shall obtain the large distance behaviour of the two- and three-point functions of 
the model. As a preliminary step in this direction we shall first obtain an explicit expression for the 
two-particle form factors. 

3.1 Form factors of order /disorder operators in three-state Potts model 

The continuum limit description of the model in which we are interested is given by the thermal 
perturbation of the above discussed CFT, i.e. 

S = 5cFT + T j d?x e{x) (6) 

and belongs to the class of integrable QFTs [8] . 

In the computation of form factors it is fundamental to establish which is the nature of the basis 
of asymptotic states of the theory, i.e. which are the particle excitations which come into play. This 
is the main difference between the high-temperature and low-temperature phase of the model. 

At T < there are three degenerate vacua and the excitations are the kinks Kjj±i, with j = 1, 2, 3, 
which interpolate between the ground state j and the ground state j + l(mod3). Form factors in this 
regime can be found in [9]. 

At r > the ground state is unique, the S3 ~ Z2 x Z3 symmetry is unbroken and its simplest 
realization is in terms of a doublet of particles A and A transforming as 

nA = u;A, nA = uj^'^A; CA = A (7) 

where Q is the Z3 generator, C the charge conjugation operator and uj = e^'^*/^. The Z3 symmetry 
allows the existence of the fusion process 

AA A. (8) 
The two particle S'-matrix for this integrable model turns out to be [10, 11] 

\A{9i)A{e2)) = u{9)\A{e2)A{9i)) 

\A{9i)A{e2)) = m\A{92)A{ei)) (9) 

where^ 6 = 61 — 92 and 

nW=t(..-^)="-^f(^ + y. (10) 
' ' sinhi(e- ^) ^ ' 

The pole present in the amplitude u{9) located at 6* = ^ corresponds to the bound state (jEJ. 

For the reasons discussed above, in the following we shall concentrate on the r > case. The 
symmetry of the model implies the presence of a doublet of spin operators a and a which form a 



^The rapidity variables 6i parameterize energy and momentum of the particles of mass m as (Pi,Pi) 
(m cosh 6i , m sinh 6i ) . 
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two-dimensional representation of S3 and transform as A and A in Q. Duality also requires the 
existence of a doublet of disorder operators /i, fi. 

The n-particle form factors of a given operator <I> are defined as 

i^at...,a„ (01, . . . , 0n) = (0| <J>(0) |yl,, (^i) . . . A^^ (9^)) (11) 

with the indices Oj referring to particle species. They satisfy a set of axioms which can be used to 
compute them explicitly [12-14]. 

In the following we shall be interested in the two-particle form factors for the a (a) and /u (fl) 
operators [9,15,16]. The spin operator satisfy (we recall that a and /u have non-zero form factors 
upon charged and neutral asymptotic states, respectively) 

FUi0) = u{e)F^A0 + 2ni) (12) 

with the residue condition due to the bound state pole of 

-iRes,^2^FI^(e) = d^F| (13) 

where F^ is the one-particle form factor and the three-particle coupling constant is given by 

i Resg^2^u{e) = {ri^f = V3. (14) 

Since the disorder operator /x and fl are non-local with respect to the spin operators a and a which 
create the particles A and A, a phase factor g^^'^*/^ enters in the form factor equations 



F^^{e) = m e-2-/3^M^(0 + 2vri) 

F^^{e) = me'-'/'F^^{e + 27ri). (15) 



Defining the vacuum expectation value of the disorder operator fi as (fi) , the residue conditions due 
to the kinematic pole follow 



-iRese=^^F^^ie) = (l - e'^^/s^ (m). (16) 

In order to work out the solutions, we introduce integral representations for both u{9) and t{9). Let 
us define 

f tQ\ sinhi(6l + i7ra) 

Ja\P) 



sinh i(0 — lira) 



= exp 

Then, we have immediately 



2 

dx sinh(l — a)x , Ox 

smh — 

X sinhx iir 



(17) 



n{e) = -f2rm, m = h/M- (18) 
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Since the solution to the functional equation Fa{6) = fa{0)Fa{6 + 2i:i) admits the integral represen- 
tation 



„ , dx sinhfl — a)x . o Utt — 6)x 

2 / ^ — — — sm'' — 

27r 



X smh X 



one finds the desired expressions for the 2-particle form factors 

i sinh 0/2 



sinhi(0 + ^) sinhi(0- ^) 2^2/3 (2vri/3) 



^/3 
^/3 



(19) 
(20) 



(21) 



The symmetry of the model requires that the remaining form factors for the conjugated operators 
are given by 



(22) 



AA 



AA^ 



AA^ 



(23) 



The relative normalizations between order and disorder form factors can be obtained by means of 
the cluster condition [9, 17] 



FIF% 



(24) 



(symmetry requires F| = F^). 

The knowledge of the first few form factors gives access to approximate expressions for correlation 
functions. Their explicit expressions for the operators a, a, jj, and Jx are listed below 



a{x)a{fS)) = -^^KQ{m\x\) + ^ I AO Flj^{2e)Flj^{-2e) KQ{2m\x\cosh.e) + . . . 

TT Jo 



(25) 



1 



(//(a;)/2(0)) = ii^y + 



de F^^{29)F^^{-2e) + F^^{29)F^^{-29) Ko{2m\x\coshe) + . . . 

1 f'°° r 1 
(M^)MO)) = {f,f + -J^ de [F'^^{2e)F'^^i-2e) + F'^^{2e)F^^{-2e)\ Ko{2m\x\coshe) + ... 

where ellipsis stand for terms of order e"^'"!^! for large \x\. We remember that symmetry requires 

(m(x)/z(o)) = 

(m(x)//(0)) = (/i(x)/2(0)). 
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3.2 Three-point correlation function 

As first noticed in [19], the pole due to the process ((SJ has a pecuhar effect on the spectral expansion 
of the three-point correlation function G^^\xi,X2,X3) = {a{xi)a{x2)o-{x3)) . We recall that (in the 
notations of high-temperature phase) it takes the following form 

G(3)(xi,X2,X3) = ^^21!£m j^o(^r^)+o(e-"0 (26) 

vr 

where it is assumed that all the angles of the triangle {xi,X2,X3) are less than 27r/3 and ry denotes 
the minimal total length of lines connecting the 3 spins and is given by 



ry = n + r2 + rs = y ^Rj^ + + ^^li + 2^^5a • (27) 

Here 5a is the area of the triangle {xi,X2,X3) and the meaning of and Rjk is clear from figure^ 
The point at which the three segments of lengths rj join is usually called "Steiner point". If one of 
the angles becomes larger than 27r/3, then the large distance asymptotic is dominated by 0{e~"^^) 
where 

p = min(i?i2 + i?23 , ^23 + R31 , R12 + ^31 ) • (28) 

The aim of this section is to give a proof of the previous result together with an explicit expression 
of the subleading term 0{e~"^P). 

The spectral expansion of the correlation function G^^'^ can be obtained inserting the resolution 
of the identity between the operators which appear in the correlator. More explicitly, we can write 
the general expansion as 



Gf'\xi.X2,X3) = 53(0|<r(ii)|fc)(fc|»(i2)|/)(i|»(i3)|0)=j;og„(j!i,i2,i3) (29) 



k,l k,l 



and hence we will use the notation [k, I] to identify a given contribution (k and I denote the number 
of particles present in the states \k) and |/), respectively). 

It is understood that all the following formula are valid in the limit of large triangles, namely 
when the condition mRi2, mi?23) fnRvi ^ 1 is fulfilled. 

3.3 [1,1] contribution 

The leading term of the spectral series is 

l-CO JO l-CO JQ 

GjJ),j(xi,X2,X3) = y ^-(O|a(xi)|i(/3))(A(/5)|a(x2)|A(0))(A(0)|a(x3)|O). 

Extracting the space-time dependence of the matrix elements the previous expression becomes (see 
figure n for conventions) 

/oo ^Q l-CO JQ 
2^ y 2^^^^^^ -e + in) exp [-im {Ru sinh(/? - ta) - R23 sinh(0 + ^7))] . (30) 

Then, taking into account that the residue on the bound state gives the condition 
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Figure 1: Triangle associated to the three-point correlation function. The orientation of the triangle 
has been chosen so that one of the three segments joining the vertices with the Steiner point lies on 
the real axis. 
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(a+ = a + 7) and using the Cauchy's theorem (see Appendix A for details), we can write Gyi i^ as 
follows 

4'V^l'^2,^3) = ^ (y - "+) ^^^^^^i^o(mry) + (31) 

/oo />00 JQ 

TT / - e + ^(a + 7)) exp [-m (i?i2 cosh/3 + R23 coshO)] . 

-00 ^'^ J-00 ^'^ 

where -ff (z) is the usual step-function and ry = J2i=i (see fig. 
3.4 [1,2] & [2, 1] contributions 

According to the general formula, the next-to-leading term of the spectral series is given by 

1 /-oo sin POO J/3 _ _ _ _ _ _ 

I I 1 '(O|a(xi)|i(/3))(i(/3)|a(x2)|i(0i)i(02))(A(02)A(^i)|cT(x3)|O) + 



1 /-oo J/O /"OO J/1 ^/D 

+ 2! y y ^^^(O|a(xi)|^(ei)yl(02))(A(^2)A(0i)k(x2)|A(/3))(A(/3)|cT(x3)|O) 

and hence, extracting the space-time dependence of the matrix elements, we obtain 

Gap (X1,X2,X3) = G^^^^^^{xi,X2,Xi,)+G^^^^-^{xi,X2,X'i) = 

1 p d/? ^61^62 



(O|a(O)|A(/3))(A(/3)|(7(O)|A(0i)A(^2))(A(^2)A(^i)k(O)|O) x 



2! 27r (27r)2 

exp [—im {R12 sinh(/3 — ia) — i?23(sinh(0i -|- i^) + sinh(6'2 -|- ^7)))] + 

21 2^ 72^ 

exp [—im (i?i2(sinh(0i — ia) + sinh(02 — io)) — R23 sinh(/3 -|- i^))] ■ (32) 



/ I -^{O\amAi9^)Ai02)){Ai92)Ai9^)\amMP)){Am^m\O) x 



It turns out that the above term gives a contribution to the correlator of the same order of magnitude 

(3) 

of G'ji !]• In order to understand this let us first notice that the matrix elements have to be rewritten 
in terms of form factors (which are the quantities that one actually computes in integrable field 
theories). Such a reduction to form factors is achieved by iterating the crossing relation [12, 13,20], 

(3) (3) 

and in the case of G^i 2] (^^^ treatment of ^] will follow the same guidelines) we have 

{A{(3)\amA{9i)A{e2)) = {0\amA{(3 + i7T)A{ei)A{92))+ (33) 
+2Tr6{(3 - 9i){0\amA{f3)) + 2^5(/3 - ^2)^/(^1 - /3)(0|a(0)|A(/3)) = 
= FIaaH^ + ^1' ^2) + 2tt{8{/3 - 9i) + 6{/3 - 92)u{9i - /3))F| . 
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(3) 

Plugging the previous expression in G^i 2] performing the integration over the delta functions, 
we obtain 

pi rco JO i-co jn jq 

g[5](xi,X2,X3) = ^£ J ^^^F-^^{P + i7r,9r,92)FM02-ei) x (34) 

exp [—im (i?i2 sinh(/3 — ia) — i?23(sinh(f?i + i^) + sinh(6'2 + *7)))] + 
+^-^ J -{^^aa{02 - Oi) exp i-im {Ru smh{9i - ia) - i?23(sinh(ei + i^) + sinh(02 + n)))] + 



1 J 7^^^^^^' " '^''P ~ '"^ ~ -R23(sinh(ei + i7) + sinh(02 + il)))] 



where in the last line we have used the fact that F'^j^{9) = tt(6')F^^(— 0). A brief inspection of such 
an expression shows that the last two terms coincide (to see this one can make the exchange 0\ <r^ 62 
in one of them), and hence we can write 

Gf^^^{x,,X2,x^)= j ^^j ^-±-^Fl^^{P + m,0,,e2)nA{e2-e,)x (35) 

exp [—im {R12 sinh(/3 — ia) — /223(sinh(^i + ij) + sinh(02 + ^7)))] + 

/OO Jfl J/3 
--L-^F^^{92 - 9i) exp [-im {Ru sinh(0i - ia) - R23{smh{9i + ij) + sinh(^2 + il)))] • 
-00 V^^j 

(3) 

Such an expression for Cj]^ 2] shows that the disconnected parts sum up giving a contribution which 
is quite similar to G^\j ■ Hence, if we pose 

G[J)^f ^^xi, X2, xs) = {F^)' 1^ ^^Ia(^2 - ^i) X (36) 

exp [—im (i?i2 sinh(0i — ia) — i?23(sinh(0i + ^7) + sinh(02 + ^7)))] 

and we use the Cauchy's theorem taking into account the residue condition on the bound state, we 
obtain (see Appendix A for details) 

G[gf ^^xi, X2, X3) = - f ) ^SllAAK.imry) + (37) 

/OO JQ ^0 
-J^^aa{02 -9i + #) exp [-m (i?23 cosh 92 + Riz cosh 9^)] . 

One can notice the change of sign in front of the function H{^z) and the difference in its argument 
with respect to Gj^^^^j. 

As anticipated, we can treat ™ exactly the same way as 2] • °^ computation 

we are left with the following expression for G^^^^^'^ 

Ggif^7xi,X2,X3) = -h[4>-'^^) ^I^lAAK,{mTy) + (38) 

/OO JQ JQ 
-J^FMOi -O2 + i(t>) exp [-m {Ri^ cosh 9^ + R12 cosh ^2)] • 
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Collecting all the previous contributions we can write down the explicit expression of the spectral 
series for the three-point function up to one-particle contributions 

(3) ^{3),disc ^{3),disc 



^27r \ 27r\ / , 27rM ^^^1^ 



/oo ^ fjAf) 
TTTTI -0 + ia+) exp [-m {Ru cosh p + R23 cosh 0)] + 
-00 l^'^j 

-O + i-fp) exp [-m (i?23 cosh /? + i?i3 cosh 9)] + 

-O + icf)) exp [-m (i?i3 cosh /? + R12 cosh 6*)] + . . . . (39) 

Some comments about the previous expression are in order. 

• The combination of H[z) functions which appears in the first term nicely realizes the condition 
stated at the beginning of the section. It is not difficult to show that only two possibilities are 
allowed: either all the angles are less than 27r/3 and such a prefactor is 1, or one of them is 
larger than 27r/3 and then the prefactor is zero. 

• The second term of ()39() shows explicitly that its asymptotic behaviour for large distances is 
given by 0{e~'^P) where 

p = min(i?i2 + i?23, ^23 + -R31, + ^3l)- 
3.5 Isosceles S>c Equilateral triangles 

Formula ()39() can be simplified by means of a wise choice of the geometry of the triangle. 

Let us begin with the case of isosceles triangles. We decide to make the choice: R = Ru = R23, 
S = i?i3 = 2R sin a, and = 7 (a+ = 2a), -0 = 1;^ = 7r/2 — a (it is easy to show that any other choice 
of the isosceles triangle will be equivalent to the previous one). Hence (|39|) will take the following 
form 

G(3)(xi,X2,X3) ~i/(|-a) ^^^^^J^'^ Ko (2mRcos (f - «)) + 

/oo jn / n\ 

^FU{9 + 2ia)Ko (^2mi?cosh-j + 

+2(F|)2 / ^^^Fl^ifi -9 + i{7r/2 - a)) exp [-mR (cosh /? + 2 sin q cosh 9)] . 
J-00 (27r)^ 

An even simpler form can be found when the case of equilateral triangles is considered. Such a choice 
corresponds to fixing a = vr/G in the previous formula (obviously one can obtain the same result 
directly from formula (|39j) ) finding 

G(3)(xi,X2,X3) ^ ^-^L_AAko (V^mR) + 

^ (2^^^^^^ + '""/^^^^ (^2mi?cosh -J . (40) 
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Figure 2: Triangular lattice. 

4 Monte-Carlo simulations 

We performed Monte-Carlo simulations of the 3-state Potts model on a triangular lattice (where each 
site has six nearest neighbors, see fig. |2) with honeycomb boundary conditions. This lattice has a 
larger symmetry than a square lattice and allowed us to study 3-point functions with three exactly 
equidistant points. A Swendsen-Wang cluster algorithm was used and 2- and 3-point functions were 
computed from improved estimators. The latter proved essential for the 3-point function, which 
could not have been measured to as large distances with conventional estimators. Simulations were 
performed at three different couplings, see table^below. In the table, a volume 60 x 180, for instance, 
refers to a hexagon with side 60. 

4.1 Determination of the mass and F| 

The long-distance expansion for the 2-point function has two free parameters, the mass (inverse 
correlation length) and the 1-particle form factor. For each lattice, the mass was extracted from the 
effective mass, defined in terms of the wall- wall correlator Gq{x) as^ 

meff(x + i)^21n^^^^ (41) 

On a triangular lattice, 

G'oix) = ^(f7,(x, y)ae{0, 0)) ~ ^ dy {a,{x, y)a,(0, 0)) (42) 

where the factor comes from the distance between points along the wall on a triangular lattice, 

see fig. 121 Note also that a displacement of ^ instead of 1 is possible on such a lattice. 

The prediction of 1- and 2-particle contributions for the propagator including finite-size effects 
for the 1-particle contribution (those from the 2-particle contribution are negligible) is 

Gl{x) = fe""*^ + g2{mx) + e''^^^"^)) (43) 

V3m ^ ^ 



where 



oo 



g2{mx) = ^ I de 







exp[— 2mj; cosh ^1 , , 

coshg 



^We denote lattice quantities which differ in normahzation from their field theory counterparts by an index 
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13 volume ^ m Fi Xmin a^max r4 



AA 



0.5900 60 X 180 4.5184(9) 0.22132(5) 0.92179(19) 7.5 25 1.2609(16) 
0.6130 120 X 360 9.093(3) 0.10997(3) 0.8233(2) 15 50 1.287(3) 
0.6231 240 X 720 18.232(7) 0.054849(19) 0.7420(3) 25 125 1.299(3) 

Table 1: Simulation parameters. 

and Fi is the 1-particle form factor F| in lattice normalization. The corresponding effective mass is 

meff(a; + i) = m + Jmcfr(x + i,m,L) (45) 

with 

orripff x + -7, m, = 2 in -. — -— ^ tt — - — (4d) 

^ii^ ^4' ' ^ l + e™(^+i/2)g2(m(x + i)) + e-'"(^-2^-i) ^ ^ 

The fit with this function is performed by subtracting 5megF(x, m, L) with some trial m from the 
effective mass data, fitting the result to a constant, and then iteratively improving m. This procedure 
is very stable because the corrections are small in the range of the fit. The fit at /3 = 0.613 is shown 
in fig. 121 It confirms the analytic prediction very well. The resulting masses and correlation lengths, 
together with the fit ranges used, can be found in tabled The error estimates were obtained with a 
jackknife analysis. 




10 20 30 40 50 60 70 80 90 



X 

Figure 3: Effective mass at /3 = 0.613 and 1-parameter fit from eq. (|45|) including the long-distance 
expansion up to 2 particles and finite-size effects, x has period 180. 

Given the mass, the 1-particle form factor can be obtained from the wall-wall correlator by fitting 
the prefactor in (|l3|l. The result at /? = 0.613 is shown in fig. |3 The long-distance expansion and 
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finite-size effects describe tlie data very well. The resulting 1-particle form factors are given in 
table n The error estimates were again obtained from a jackknife analysis which also included the 
determination of the mass described above. 




3.54 - 



10 20 30 40 50 60 70 80 90 

X 

Figure 4: Wall-wall correlator at /3 = 0.613 and 2-particle prediction with finite-size effects. Only 
the prefactor has been fitted. 

The long-distance expansion of the 2-point function is completely determined by m and Fi. 
Figure [5] shows the 2-point function divided by the 1-particle term of the long-distance expansion 
and rescaled by m and F^. The data from the three different lattices collaps nicely, and the small 
deviations from the 1-particle term are well described by the 2-particle term. 



4.2 Scaling violations 

Comparing the results from the different lattices, we may estimate the size of the scaling violations 
which we must expect in our results. 
Setting 

we expect m ~ 

t^/e and F/, ~ t^/^. Fi gures ini and [7| show the ratios m/t^/^ and Fn/t^/^ as a function 
of t. They vary by about 3%, a reasonable size for scaling violations at these correlation lengths. 
The scaling violations in the mass are expected to be of the form 

m{t)/t''^ = mc{l + am,At^ + am,i* + •••)• (48) 
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Figure 5: Two-point function, rescaled using m and Fi from fits to the wall-wall correlator and 
divided by 1-particle term of long-distance expansion. 
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Figure 6: Correction to scaling for m. 
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Figure 7: Correction to scaling for F^. 

Taking into account that the subleading thermal operator e' corresponds to the conformal operator 
(3, 1; 3, 1), the leading correction-to-scaling exponent is A = — 2-x' — |- This value has been first 
determined in [21] and later confirmed in [22-24], see [25]. We find am, A cLm,i- A fit with only one 
term and variable exponent p yields p = 1.10 ... . These two facts indicate that am, A may vanish. A 
fit with only the linear term has = 7.7, so we perform a fit with the next two terms, 

mc(l + am,it + am,2A't^^) ■ (49) 

Since it has as many parameters as data points, we use the fit with only the linear term for an error 
estimate. The result is 

rric = 2.112(3) . (50) 

The statistical error from the fit with only a linear term is much smaller. The fits are shown in fig.lHl 
For Fi, a fit to 

Fi{t)/t^/^ = Fe(l + aF,At^ + aF,it + ...) (51) 

yields a linear term which is much smaller than the leading correction. We use a fit with only the 
leading correction for an error estimate and obtain 

Fc = 1.189(3) . (52) 

Again the statistical errors from the latter fit are smaller than the quoted error. Both fits are shown 
in fig. El 
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4.3 3-point function 

The long-distance expansion of the 3-point function is completely determined by m and F^. In 
particular, G^p /Ff as a function of mxi is independent of any parameters. The measured 3-point 
functions, rescaled in this way, are shown in fig. |H1 Also shown is the long-distance expansion of eq. 
(|in|l which we split into "Y-type" and "A-type" contributions (the second term depends on sums of 
two spin-spin-distances) , 



Gf\xi,X2,Xz) = FfTj^ji\gY{xi,X2,X'i) + gK{xi,X2,Xn)\ 



with 



C/y(2;i,X2,X3) 
S'a(2;i,X2,X3) 



-Ko{mrY) 
vr 

3 
2T2 



(53) 

(54) 
(55) 



where R = \xi — X2\ = 1^2 — X3I = |x3 — xi| is the distance between two spins and ry = ^/SR the 
Y-length. For the equilateral geometry, the Y-type term is leading. It is shown separately in the 
figure. The subleading A-type term is sizable up to large distances. This is not unexpected as the 
sum of two sides, 2R, is only 15% larger than the Y length \/3R. 
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Figure 8: Three-point function, rescaled using m and F^ from fits to the wall-wall correlator. 



The deviation from the leading Y-type term is shown in fig. |3 Scaling violations are clearly 
visible and turn out to be of the same order as those observed in the 2-point function, i.e. about 
3%. These scaling violations have to be expected (even after rescaling with m and Fi) since the 
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lattice model at a given temperature t is effectively described by a quantum field theory which 
differs from the scaling theory © by irrelevant operators. The form factors therefore differ from the 
values computed exactly in the integrable scaling theory. The two terms above contain the three- 
particle coupling constant and the 2-particle form factor F^^{6) which is also proportional to 
r^^. Lacking an understanding of the rapidity-dependence of the scaling violations, we proceed by 
ignoring all scale-dependence in -F^^ other than that induced by T^^, i.e. we replace the prefactor 
in H53() by a function of t. This works very well, as can be seen from the figure. The remaining 
deviations may well be due to higher-order terms in the long-distance expansion. The resulting values 
for the three-particle coupling are given in table ^ The quoted errors were taken from the 3-point 

(3) 

function at distances where /{dY + 5a) approaches a constant. They are therefore only rough 
estimates. The values of r^^(t) nicely approach the value 3^/^ = 1.316 ... in the scaling theory. The 
corresponding curve is also shown in the figure. Note that this curve does not depend on any fits, 
but is completely determined by the values of m and extracted from the wall-wall correlator. 
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Figure 9: Same as fig. El but divided by the leading (Y-type) contribution. Note the larger range. 

5 Short distance expansion 
5.1 2-point function 

In order to have good control on the (era) correlation function it is worth to compare the Monte- 
Carlo data in the short distance regime with the corresponding perturbative expansion obtained in 
the framework of the so-called Conformal Perturbation Theory. Since such a perturbative expansion 
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is expected to be valid in the region m\x\ <^ 1, it is a complementary tool with respect to the form 
factor expansion discussed before. 

Let us recall the main results about Conformal Perturbation Theory in the special case of the 
correlation function {a{x)a{0)). Following the standard literature on the subject [26]- [29] we may 
write such a correlator as 

G,{x) = (a(x)a(O)) = ^cit^(x;r)([<^^]) (56) 

p 

here the sum over p ranges over all the conformal families allowed by the Operator Product Expansion 
of a and a. The Wilson coefficients C^J'^^ can be calculated perturbatively in the coupling constant 
r. Their Taylor expansion (w.r.t. r) is: 

*k^;r) = 5:^5.^cS''(x;0). (57) 

k 

It is possible to show that the derivatives of the Wilson coefficients which appear in the previous 
expression can be written in terms of multiple integrals of the conformal correlators 

d^C^t'\x; 0) = (-1)*^ J d^zi . . . d^Zk {a{x)a{0)e{zi) . . . e{zk)[^p]{^)U (58) 

where the prime on the integral implies a suitable treatment of the IR divergences and we used a 
shorthand notation to denote the follwoing limit: 

/ f \ r w r (Q"!^!) e(w))rft , . 

{(^{xi) e(oo)),ft = hm , . . • (59) 

w-^oo (e(u;)e(0))eft 

For all the details we address the interested reader to the original literature [27] . 

Another important piece of information which enters the perturbative expansion of the correlator 
is represented by the Vacuum Expectation Values ([0p]). Such VEVs cannot be calculated in the 

framework of the perturbation theory, and being non-perturbative objects they have to be obtained 
by other methods. By using the powerful techniques of Integrable QFTs, they were computed in a 
series of papers [30]- [32] for a wide class of theories, including various integrable perturbations of 
the Minimal Models. 

The leading term of the perturbative expansion is given by the 2-point conformal correlator which 
corresponds to the choice A; = and 0p = I 

cUx;0) = j-^ (60) 

which follows from the choice of the conformal normalization C^^ = 1. 
The first few sub-leading terms are 

G,{x) = (a(x)a(O)) =C(x;r) -FC^^(x;T)(e) + ... (61) 

where 

CI^{x;t) = Cl^{x;0)+TdrCl^ix;0) + ... 

C^(x;r) = C^^(x;0) + ..., (62) 
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give corrections up to r. The explicit expression of the various contributions is, together with (|6flj) . 



drCl^{x;f)) 



sm ■ 



r(-i/5)r(3/5) 



r(2/5) 



where the Wilson coefficient 



1 /sinf r^(3/5) 



CI 



2\ sinf r(2/5)r(4/5) 



0.546178... 



(63) 



(64) 



can be found combining the results of [33] and [34] (see also [35]). The integral which appears in 
drC\^{x] 0) is well known. It is a particular case of the following one 



smvra 



r(-a-6 - l)r(6 + l) 



r(-a) 



(65) 



-8.97743 . 



Its value is: y_2 2 

5 ' 5 

At last, let us discuss the other (non-perturbative) quantities we need: the VEV of the perturbing 
operator e(x) and the relation between the coupling constant and the mass of the fundamental 
particle. The latter is known, and can be extracted from [36] 



r = Km^/^, fi; = 0.164303, 



(66) 



The VEV (e) can be easily computed starting from the knowledge of the vacuum energy density [37] 



V3 



-m 



which is related to the VEV of the perturbing operator by means of the relation 

(e) = drSo = A, r2/3 = -9.761465 . . . r^/^ = -2.92827 . . . m^/^ 

Collecting the above ingredients, the perturbative series can be cast in the following form 

1 



(67) 



(68) 



G^x) 



|x|4/15 
1 

1 



\X 



|4/15 



5 ' 5 



where we set u = t\x\^^^, and r = m\x\, respectively. 
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5.2 3-point function 

The perturbative approach of the previous section can be generahzed to multipoint correlators as 
well. Following the guidelines of [29] one is able to write down a perturbative IR safe expansion 



p 



where the structure functions c'^^^^ {x; r) are a generalization of those appearing in the formula for 
the two-point function^. Their Taylor expansion gives 

OC £ 

c't:U = 0) + E ^rc'tt,. o) (70) 

(.=1 



where we have 



and 



4t'.<^„fe0) = (<^l(xi)..>n(Xn)[</>p](oo)),ft, (71) 



^r4tl<^„fe0) = {-lY j d^zi...dht{Mxi)...M^n)e{zi)...e{ze)[cPp]{^))^,, (72) 

where the operator e{x) which appears in the previous expression is the perturbing operator conju- 
gated to the coupling constant r. 

In the present section we are interested in the 3-point function G^^^ = (o"icr2cr3), whose perturba- 
tive expression can be written according to the previous considerations as 

G^^\x) = {a{xi)a{x2)<j{x3)) = 

= Claa{^i^X2,X3;T) +C^^^{xi,X2,X3;t) {e) + ... ; (73) 

up to first order in r one has 

(a:i,X2,a:3;0) + t drCl^^{xi, X2, X3;0) + ... 

Caaai^l^X'2,X3;T) = C^^^(xi, X2, X3; 0) + T 5T-C^^^(a;i, X2, X3; 0) + . . . 

The explicit expression of the zero-th order contributions can be derived quite easily. We have 



Claa{xi,X2,X3;0) = {a{xi)a{x2)a{x3)). 



(74) 



\XI - X2|2/15|X2 - X3|2/15|a;i - X3\^/^^ 

where the structure constant C^^ can be found in the literature [33-35] 



_ 3sin27r/5 r(5/6)r^(3/5)r^(l/3) _ 

" V 7rsin7r/5 r(2/5)r(4/5)r2(2/3)r2(l/6) " ' ' ■ ' ^^^^ 



*We used the notation x = {x\, . . . , 
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The other zero-th order term is given by 

Caaa{xi^X2,X3;0) = {a{xi)a{x2)(T{x3)e{oo)) (76) 

whose main ingredient is the conformal fom' point correlator (aaae) which was computed in [41]. We 
consistently fixed its normalization constant, which gives 

{a{x,)a{x2)a{xMx,)U = C^^CU 'g/ls Ivi^ - y)\''^'' {\fi{y)\' + K-'\h{y)\'} (77) 

|a;i4a;24a:;34r/-^^ 



where 



and 



/i(y)=y~'/'2Fi(l/5,4/5;2/5;y), ^(y) = 2^1(4/5, 7/5; 8/5; y) (78) 



^ 9r(l/5)r(3/5)^ ^ XUX23 _ ^ X13X24. , . 

"4r(4/5)r(2/5)3' ^"X34X21' ^"X34X12 ^ ^ 

Then, a simple computation gives 

C^,,(xi,x2,x3;0) = Cf,Q^|xi2Xi3X23p/^^|y(l-y)|^^/^^ {\hiyr + K-'\f2{y)\'} (80) 

where 

- a:23 . ~ a;i3 

y = — , 1 - y = — (81) 

a;2i Xi2 

The first order terms are as follows 

drCl^„{xi,X2,X3;0) = - j (f'z{a{xi)(j{x2)(T{x3)e{z)),i, 

drCl^^{xi,X2,X3;Qi) = - (fz{a{xi)a{x2)a{x3)e{w)e{oo)),ft (82) 



but we did not manage to calculate them analytically for a generic triangle. We expect that some 
simplifications will occur when a given geometry is chosen, for example in the case of equilateral 
triangles. 



6 Comparison with Monte- Carlo data at short distance 
6.1 2-point function 

The short distance expansion calculated in section [S] can be compared with the Monte-Carlo data in 
the region of short distances, i.e. m\x\ <C 1. 

For this purpose, it is first necessary to rewrite such an expansion in the dimensionless variable 

r = m\x\ 

G^r) = m~'/''{a{x)am = ^ (l + 9ir'^' + 92r'/' + gsr' + g^'^' + . . . ) (83) 
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where, from the previous section we have 

c/i =C^^A,k2/3 = -1.59936..., 52 = 3^„2_2C^s'^ = 0-805622..., (84) 

and and (74 are unknown constants which embody the contributions given by higher orders in 
perturbation theory. 

We chose to use the sample of Monte-Carlo data with the largest correlation length, i.e. ^ = 
18.232(7). This is motivated by the following requirements: 

• The perturbative calculation is supposed to hold for r <^ 1; 

• The region of very short distances {\x\ < 5, i.e. r = < 0.27) have to be avoided because of 
the so-called lattice artifacts, which are non-universal corrections induced by the lattice when 
the distance is comparable with the lattice spacing (see [38] for an analysis of this point in the 
context of the Ising model in magnetic field). 

Such requirements imply the presence of a limited "window" in which the comparison between the 
data and the perturbative expansion can be done. Hence, the choice we made is motivated by the 
necessity of maximizing the number of data which fall in such a window in order to have a reliable 
sample for the fitting procedure. 

We used the following fitting function 

Giir) = m-'/^'{ae{x)am) = ^ + 9ir'/' + 92^^/^ + gsr' + g^r''/"^) . (85) 

where the constants gi and 52 are known analytically and can be fixed exactly, while the constants 
53) 94: which parametrize higher order contributions in the perturbative expansions are left as free 
parameters and are fixed by the fit. The overall constant which takes into account the different 
normalization of a on the lattice and in the continuum (the normalization of the latter has been 
fixed in H83|) and is the usual conformal normalization) is also considered as a free parameter for the 
fitting procedure. 

It turns out that lattice effects are larger than the statistical errors up to distances r ~ 0.5 once 
correlations between the data points are taken into account. This makes it impossible to use the 
goodness of fit as a criterion for a successful description of the data and for error estimates. Instead, 
we vary the fit range in order to estimate the effect of lattice artifacts, and include higher-order 
terms in order to estimate their effect. For our best estimate of the parameters, we use the shortest 
distance (r = 0.3) where lattice artifacts are smaller than (uncorrelated) statistical errors as a lower 
bound, and the largest distance where the fit still works well (r = 0.8) as an upper bound. We then 
shift the range down (to [0.14,0.5]) or up (to [0.4, 1.0]) to get a rough estimate of the uncertainties. 
We also performed a fit with terms up to g^. We included two more orders because the exponents 
of the next two terms are close (16/5 and 18/5), so cancelations have to be expected. The upper 
bound has been increased to the maximal value where the fit still works. The resulting parameters 
are shown in table |21 The largest uncertainties come from the lower bound, i.e. from lattice artifacts. 
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parameters 


3^min 






53 


54 


55 56 




0.3 


0.8 


0.6576 


0.1674 


-0.1297 




Na,g3-,9A 


0.14 


0.5 


0.6580 


0.1492 


-0.1069 






0.4 


1.0 


0.6578 


0.1641 


-0.1261 




^0-, 53, 54, 55, 56 


0.3 


1.5 


0.6574 


0.1687 


-0.1258 


-0.0386 0.0242 



Table 2: Results of fits to 2-point function at (3 = 0.6231. 
We therefore estimate the parameters as 

= 0.6576(4) (86) 
<73 = 0.167(18) (87) 
= -0.130(23) . (88) 

It must be stressed that the values of the constants gs, g^ could be affected by systematic errors 
larger than the quoted errors due to possible cancellations between the two contributions. The 
quoted uncertainties also do not include the effect of scaling violations as they were obtained from a 
single correlation length. The lattice artifacts on the coarser lattices are too strong to allow for an 
estimate of the scaling violations in 173 and g^. It would be very interesting to fix this uncertainty 
by a direct calculation of these coefficient, which could be performed in the framework of Conformal 
Perturbation Theory but requires techniques more sophisticated than those discussed in this paper. 

We can get a rough estimate for the scaling violations in A^^-, by fitting it to the data on the coarser 
lattices while keeping (73 and g^ fixed. Allowing for deviations due to the large lattice artifacts 
observed at small r, this works quite well. The results are N^j = 0.6725(8) for /3 = 0.613 and 
No- = 0.6999(9) for /3 = 0.59. The uncertainties include only the effect of varying the fit range. We 
have no way of estimating the bias caused by keeping and (74 fixed, let alone the effect of higher 
orders. The two fits as well as the above fit to the data at (3 = 0.6231 are shown in fig. 1101 The 
rescaled correlation function G'o-(r) has been divided by the leading short-distance term r~^/^^ for 
this figure. Since the series has fractional exponents, we use r^/^ for the abscissa. 

The scaling violations in Nu are expected to be of the form 

N^{t) = N^{1 + aNAt^ + ■ ■ ■) (89) 

with the correction-to-scaling exponent A = 2/3 that already appeared in the discussion of the mass 
and the form factor in sect. 14.21 Remarkably, our crude estimates of Nu follow this prediction almost 
without deviation. Figure ^2 shows a fit with A = 2/3 as well as a fit with two correction-to- 
scaling terms with exponents 2/3 and 1, which deviates only slightly. From the former, we obtain a 
continuum extrapolation of 

= 0.6376(6) (90) 

with the error taken as the difference of the two fits. Because of the unknown systematic errors in 
the individual N^, this error is probably strongly underestimated. 
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Figure 10: 2-point function divided by leading short-distance term versus r^/^. Fit up to to 
(3 = 0.6231 with range [0.3, 0.8] (red) and series up to with only N^^ fitted to /? = 0.613 (green) 
and /3 = 0.59 (blue) with ranges [0.4,0.8] and [0.6,1.0], resp. Theoretically known part of series up 
to g2 with prefactor from continuum extrapolation of fig. ^2 (black) . 

The short-distance expansion with this prefactor, and only including the analytically known terms 
up to g2 is shown in fig. 1101 The discrepancy between the data and this curve has two origins: scaling 
violations in the lattice data and omission of higher-order terms in the short-distance expansion. In 
the continuum limit, and with higher-order terms included, the two are expected to meet. The figure 
matches this expectation nicely. 

The constant Nu plays the role of the normalization of the magnetization operator cr^ on the 
lattice, and for this reason it is not calculable in the framework of field theory. It is a non-universal 
quantity which depends, among other things, upon the geometry of the lattice. In the case of the 
Ising model on a square lattice the analogous quantity can be evaluated exactly thanks to the fact 
that the aa correlator can be evaluated analytically on the lattice at the critical point. By comparing 
this result with the continuum-limit expression the normalization constant is then easily obtained 
(see for instance [40]). Unfortunately a similar exact lattice result does not exist for the three states 
Potts model on a triangular lattice. 

6.2 Relation of short- and long-distance expansions 

Given the normalization of the short-distance expansion N(j, the amplitude of the 1-particle form 
factor can in principle be calculated in field theory using the cluster condition (j2l|) . The explicit 
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Figure 11: Scale dependence of N^r- 



form of F^^{9) gives 



0.968240 . . . 



(91) 



The required amplitude of (//) in field theory, however, has not been computed so far. We can give 
an estimate of the latter using the measured values of A^^ and Fi. Since ai ~ \nv^a, the lattice and 
field theory form factors are related as = ^/TV^F^. 

In order to eliminate corrections to scaling, we extrapolate to the continuum limit. The scale 
dependence on the lattice was given in sect. 14.21 in terms of t, 

Fe ~ F,t^/^ (92) 
m ~ ruct^/^ (93) 

while in field theory it is usually given in terms of the coefficient r of the energy operator, 

~ Aft^^^ (94) 



(95) 



The ratio F^/m?^^^ is scale- independent. The corresponding amplitude ratio can thus be identified 
on the lattice and in field theory, 

Fjy^ Af 



2/15 



.2/15 • 



The measured amplitudes (ISU]) and together with (|Un|l yield the ratio 



2/15 

rric 



1.348(4) . 



(96) 



(97) 
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Using the known value Am = n ^/^ = 4.50431 . . . from (|66|) . we get 

Af = 1.647(5) . (98) 
The amphtude of the disorder parameter is defined via 

{^^) ~ A^rV9 . (99) 
The cluster condition (|91j) then gives A^ = Apj^fM and thus 

A^, = 1.674(5) . (100) 

All errors are probably underestimated, mainly because of the uncertainties in A'^^ as discussed in 
the previous section. 

6.3 3-point correlator: Equilateral triangles 

Let us write explicitly the perturbative expansion of the 3-point correlator in the case of the equilateral 
triangle. It will be useful in the perspective of a comparison between Monte-Carlo data and the form 
factor expansion worked out previously. 

In such a particular case some simplifications occur, i.e. \xij\ = R where R is the side of the 
triangle. Then, with the choice (obviously the final results will be independent of such a choice) 

xi = 0, X2 = i^e*^/^ X3 = i?e~^^/^ (101) 

we have 

y = i e^We^ l-y = e^^^/^ (-102) 

and finally 

C(0,i^e-/^i^e— /^0) = C7^,Q^ {|/i(ie— /6)|2 + ir-i|/2(ie-*-/6)|2}i?2/5 = Q^^i?2/5(io3) 
where 

= 0.788825 .... (104) 
Taking into account all the other pieces we have 

Such an expression is very similar to the perturbative expansion of (era), in particular one can write 
it down using the dimensionless variable u = tR^/^. As a consequence it is possible to guess the 
functional form of the higher order terms relying on dimensional analysis only. Hence we can write, 
in terms of r = mR, 

+ C4 + C5 

ri2/5 + ...) (106) 
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where 

ci = = 1.09236 . . . , C2 = C^^^ k^/^ = -2.29795 . . . (107) 

and C3, . . . are unknown (but in principle calculable) constants. It is also useful to define the scale 
invariant form of the correlator 

G(3)(r) =m-2/5(a(xi)a(x2)a(x3)) = ^ (ci + C2 + C3 r''/^ + C4 + C5 r^"/^ + . . . ) (108) 

and its corresponding expression on the lattice 

3/2 

of (r) = ^ (ci + C2 r^/s + C3 r^/^ + C4 + C5 r^^/s + . . . ) (109) 

where No- is (the square of) the normalization of the lattice magnetization operator. N^j has been 
extracted from the aa correlator in sect. 16. 11 Since the scaling violations in the 3-point function are 
expected to be different from those in the 2-point function, we do not use the values of N^j fitted to 
the individual lattices, but rather the continuum extrapolation N^. The series up to the known C2 
with this prefactor is shown in fig. I12L Note that the abscissa is the distance between two spins (in 
physical units and to the power 2/5), not the Y-length mry = \/^r which appears in the long-distance 
expansion. The figure shows that, while higher-order terms are clearly important in the range where 
we have data, the leading term in the series is compatible with the latter. 

A reasonable description of the data on the finest lattice can be obtained with the series up to 
C5 with No- and C3, . . . ,05 as fit parameters. Like in the case of the 2-point function, we use a fit 
range where both lattice artifacts are small compared to uncorrelated statistical errors (r > 0.3) 
and the fit still works well (r < 0.8). The fit is shown in the fig. 1121 Contrary to the case of the 
2-point function, the fit does not stay near the data beyond the fit range, and, not surprisingly, the 
parameters depend quite strongly on the latter, see tableOl Since the fit already includes three terms, 
we do not attempt to estimate the effect of higher-order terms by including even more terms. The 
following uncertainties therefore contain neither the effects of these nor those of scaling violations: 



= 0.6376(36) (110) 

C3 = 1.24(10) (111) 

C4 = 0.44(45) (112) 

C5 = -0.33(37) . (113) 



Still, it is quite remarkable that the fit yields a value of N^j almost identical to the one extracted 
from the 2-point function (the near coincidence of the curves at r = in the figure is not enforced!) 

Also contrary to the case of the 2-point function, the scaling violations on the coarser lattices 
cannot be described by just changing the prefactor N^j. We do not attempt fits with more parameters 
as there is not a window where both lattice artifacts are small and higher-order terms can be neglected. 

7 Implications for the three- quark potential 

In these last years much interest has been attracted by the study of the three-quark potential in 
Lattice Gauge Theories (LGT). Besides the obvious phenomenological interest of the problem, the 
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^min ^max C3 C4 C5 

0.3 0.8 0.6376 1.240 0.441 -0.331 
0.15 0.5 0.6340 1.343 -0.007 0.042 
0.4 1.0 0.6365 1.255 0.393 -0.298 

Table 3: Results of fits to 3-point function at /3 = 0.6231. 




.1 .2 .3 .4 .5 .6 .7.8.91 1.5 2 

r = m R 



Figure 12: 3-point function divided by leading short-distance term, theoretically known part of series 
up to C2 with prefactor from continuum extrapolation of from 2-point function, and fit up to C5 
to the data at /? = 0.6231 in the range [0.3,0.8]. 

three-quark potential is also a perfect tool for testing our understanding of the flux tube model of 
confinement and of its theoretical description in terms of effective string models. These models have 
been elaborated in the past years looking at the quark-antiquark potential and their extension to 
baryonic states is highly non trivial. Thanks to the improvement in lattice simulations (a summary 
of numerical results can be found in [44]- [48]), the qualitative behaviour of the three-quark potential 
is now rather well understood (for a recent review see [42]). 

• For large interquark distances the three-quark potential is well described by the so called Y 
law which assumes a flux tube configuration composed by three strings which originate from 
the three quarks and join in the Steiner point which has the property of minimizing the overall 
length of the three strings. This picture is also in agreement with what one would naively find 
using standard strong coupling expansion (notice however that due to the roughening transtion 
this is only a qualitative indication, and cannot be advocated as a "proof" of the Y law). An 
interesting consequence of this scenario is that one can use an effective string approach to model 
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the behaviour of these flux tubes and hence a "Liischer hke" correction to the potential should 
be expected. This correction was evaluated in [46] and succesfully compared with simulations 
of the 3d Z3 gauge model in [42] . 

• At shorter distances a smooth crossover toward the so called A law is observed. According to 
the A law the three-quark potential is well approximated by the sum of the three two-quark 
interactions. More precisely the A law assumes that the three-quark correlator (let us call it 
G3{xi,X2,xs) where Xj denotes the position of the j*'^ quark) is related to the quark-antiquark 
correlator Go( follows: 

G3(2;i,X2,X3) ~ \/G2{xi,X2) G2{x2,X-i) G2{xi,X^) (114) 

thus leading to a potential which increases linearly with the sum of the three interquark dis- 
tances. The scale where the transition between these two behaviours seems to occur, according 
to the most recent simulation is around 0.8 fm. 

To improve our understanding of the baryonic states it would be important now to have some 
quantitative insight in the above described picture as well as to have some theoretical argument to 
explain why instead of having a single shape stable for all the interquark distances a A — > y crossover 
occurs. Moreover, since the crossover region happens to occur exactly in the range of distances which 
is interesting from a phenomenological point of view, it would be important to have some kind of 
theoretical description of this crossover with which to compare the numerical data. 

In this respect the present study of the three-point function in the 2d Z3 Potts model is a perfect 
laboratory to address this problem. Besides the obvious similarity of the two settings it is also 
possible to find a direct relation, since by dimensional reduction the behaviour at high temperature 
of the three-quark correlator for a 5'C/(3) or a Z3 gauge model in (2-1-1) dimensions is mapped into 
the behvaiour of the 2d Z3 Potts three-point function (in analogy to what happens for the quark- 
antiquark potential which is mappend onto the {era) correlator). 

This mapping becomes exact in the vicinity of the deconfinement transition thanks to the fact 
that both the deconfinement transition in the SU(3) LGT in (2-1-1) dimensions and the magnetization 
transition of the three-state Potts model in two dimensions are of second order. Then, according 
to the Svetitsky-Yaffe conjecture [49], the two critical points must belong to the same universality 
class and we can use the three-state Potts model as an effective theory for the SU(3) LGT. In this 
effective description the Polyakov loops of the LGT are mapped onto the spins of the Potts model, the 
confining phase of the LGT into the high temperature phase of the spin model while the combination 
fj/T = (jNt (where T is the finite temperature of the LGT model which is equivalent to the inverse 
of Nt- the lattice size in the timelike, compactified, direction) is mapped into the mass scale of the 
spin model (i.e. the inverse of the correlation length) and sets the scale of the deviations from the 
critical behaviour. It is exactly this scale which separates the A law behaviour from the Y one. 

At the critical point (i.e. when the correlation length goes to infinity) the A law is exact. In 
fact looking at ea.()60|) and (|74() for the conformal two- and three-point correlators we see that the 
relation (|114j) is fulfilled exactly. On the other side, when the distances among the spins in the 
correlator are much larger than the correlation length, simple strong coupling arguments suggest 



30 



that the dominating configuration in the partition function must be the one which minimizes the 
distances among the spins and the Y law appears. In field theory this behaviour for large separations 
among the spins is a direct consequence of the particle fusion process (jSJ- Our analysis allows to 
follow in a rigorous way the crossover between the two limiting behaviours. 

An interesting and non trivial application of our results is that they can give some insight on 
the high-temperature regime of the string corrections (the opposite of the one studied in [42]). This 
regime is reached when the interquark distances are much larger than the size of the lattice in the 
time direction and thus coincides (following the Svetitsky YafFe mapping discussed above) with the 
large distance limit of the three-point correlator in the Potts model. A remarkable and not trivial 
feature of our result in this limit is that (when the Steiner point lies inside the triangle formed by the 
three spins) the dominating exponential behaviour is dressed by a pre-exponential factor (ry)~^/^ 
which is encoded in the Kq modified Bessel function which appears in ea.()39|). This is exactly the 
same behaviour of the two-point function in this limit and, as for the quark-antiquark case [50,51], it 
indicates that in this limit the effective string corrections give a contribution proportional to log(ry) 
(analogous of the log(r) term, with r the interquark distance in the quark-antiquark case [50-52]). 
This is a rather non trivial result, which severely constraints the possible effective string pictures for 
the three quark potential and was indeed observed, in the free effective string limit in the case of the 
Z3 lattice gauge theory in three dimensions [53]. 
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Appendix A 

In order to show explicitly how to obtain equation (|31j) . let us start from (|3n|) and perform the change 
of variables (3 ^ (3 + ia — iTT/2, 9 ^ 9 — i'j + ivr/2, so that we obtain 



(F 



'OO—ia+iiT /2 J — 00+17— i7r/2 

Another change of variables, 9± = P ziz 9, gives 



- 6* + «(a + 7)) exp [-m {R12 cosh /3 + R23 cosh 6*)] . 



27r 



27r 



X 



00— 2a-|-4-27r 

xF4^(0_ + ^«+) exp 
where a± = a it 7. Now, we can integrate the function 



m { R12 cosh 



+ R23 cosh ■ 



= FMO- + exp 



m I i?i2 cosh 



+ R23 cosh ■ 



on the complex 0_-plane along the contour Cs depicted in figure El 

£ f{9_) d9_ = Us + Ds + Rs + Ls = -2m H - a+j Res,_2^_.„^/(0_) 



(115) 



(116) 



'Cs 

where H[z) is the usual step-function and 

cS—iaj^+i-K 



Rs 



Us 



S—iaj^-\-i7r 
S—ioiJ^-\-i7r 



d9- f{9^), 



d9. f{9^) 



Ds 



d9- f{9^ 



d9. f{9.). 



-S J S~ia++iTT 

The contribution of the pole is computed using the residue on the bound state 

- i Res,_2p_.^^FI^(e_ + ia+) = T^^F^. 
Then, taking the limit S" — > 00, both Us and Ds vanish, and we are left with the following expression 



°° d/3 /""^ d^ 



-00 2vr 27r 



Paa{(^ -9 + i{a + 7)) exp [-m {R12 cosh/3 + R23 cosh 6*)] 
°° d9. 



^ AA 



27r 



X exp 



-m ( R12 cosh ( ^ + y - j + i?23 cosh ( ~ y + ^7 



The desired result can be obtained in the following way. First, the expansion of the argument of the 
exponential function gives 

f9+ i-K \ ^ [9+ in 
TZ = R12 cosh ( + Y ~ ) + -^23 cosh I— ^ + ^7 



+i 



.IT vr \ / vr vr 

R12 ( cos — cos a + sin — sin a 1 + R23 ( cos — cos 7 + sin — sin 7 

o o / \ o o 

vr vr \ ( vr 

iti2 ( sin — cos Q — cos — sin a 1 — it23 ( sin — cos 7 — cos — sin 7 



cosh — + 

sinh — , 
2 ' 
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a+ < 27r/3 



\ 


in — ia^ 
\ 




~~~~~ 27ri/3 — ia+ 



a+ > 27r/3 



Figure 13: Contour of integration in the 9- complex plane. 



then taking into account the following geometric identities (figure ^) 



we finally have 



sma sm7 
sm g sm -g 

r2 = Ki2 cos a — ri cos — = K23 cos 7 — cos — 

3 3 



= (r-i + r2 + ra) cosh ^ = ry cosh ^ 



and the term proportional to sinh ^ vanishes. Hence, by means of these simplifications we are able 
to obtain formula (jHT|) . Furthermore, with the same procedure it is possible to compute both eq. (|37|l 
and eq.(jnH|l. 



Appendix B 

A non-trivial check about the correlation functions written in section lTT] (and as a consequence, about 
the form factors) is given by the computation of universal ratios. In particular, we are interested 
in those ratios which involve the amplitudes of the high-temperature (r_|_) and "longitudinal" (r_) 
and "transverse" (Fx) low-temperature magnetic susceptibilities^. These ratios have been calculated 
in the continuum by form factors techniques in the kink basis of the low-temperature phase of the 
model with the result [9, 54] 

^ ~ 13.85, ^ ~ 0.327. (118) 

^It is worth recalling that the scaling behaviour of the susceptibilities is given by 

x+^r+T"-', x-~r_|rr^ XT-rTlrp^; (117) 

where 7 = 2 = f . 
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These values have been confirmed by lattice computations in [25] and [18]. The most accurate lattice 
estimates come from the latter paper and read 13.83(9) and 0.325(2), respectively. 

Once duality is used to relate correlators in the two phases, the QFT resuls cannot depend on 
the regime used to compute the form factors. Hence, the results (|118|1 must be reproduced in terms 
of the high-temperature form factors of the order and disorder operators we have used in this paper. 
Let us recall the definitions of the susceptibilities in terms of the correlators: 



Here 



X+ = J d X (<Ta(x)(Ta(0))r>r, 

Xt = J d^x {Op\ao,ix)ao,{0)\Op)T<n, a ^ (5 . (119) 

f^a(^) = '^sW.a - ^, s(x) = 0,1,2 (120) 

and |0^), /3 = 0, 1,2 are the vacua of the model in the low-temperature phase. The spin variables a 
and a = cr* we used in this paper read 

C7(x) = e^^"^')^, s(j;) = 0,l,2 (121) 

so that 

a^{x) = \ \a{x)e'°'^ + a{x)e'^^] . (122) 
The correlators entering ()119() can then be expressed in terms of those computed in this paper as 

(a,(a;)a«(0))r>T. = ^ (a(x)a(0))T>r. (123) 
{0,3Wa{x)aam0f^)T<n = ^ {Kx)m + COS /i(x)^(0))T>r. . (124) 

In the last equation duality has been used to trade a correlators at low-temperature for n correlators 
at high-temperature. We can now use the form factors to compute the susceptibilities up to 
two-particles 

X+ ^ N ^ 0.44674 
X- ^ ^ 0.0312375 

Xt 0.0102246 (125) 

with M given in (|91j) . Taking the ratios reproduces the results (|118|) . as expected. 
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